Spontaneous deformation of the Fermi surface due to strong correlation in the 

two-dimensional t-J model 
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Fermi surface of the two-dimensional t-J model is stud- 
ied using the variational Monte Carlo method. We study 
the Gutzwiller projected d-wave superconducting state with 
an additional variational parameter t' v corresponding to the 
next-nearest neighbor hopping term. It is found that the finite 
t' v < gives the lowest variational energy in the wide range 
of hole-doping rates. The obtained momentum distribution 
function shows that the Fermi surface deforms spontaneously. 
It is also shown that the van Hove singularity is always located 
very close to the Fermi energy. Using the Gutzwiller approxi- 
mation, we show that this spontaneous deformation is due to 
the Gutzwiller projection operator or the strong correlation. 
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The effect of strong correlation is one of the most 
important issues for understanding the high-T c super- 
conductivity (SC). Among various anomalous electronic 
properties, the experiments of angle resolved photoemis- 
sion spectroscopy (ARPES) have revealed that a flat 
band around (tt, 0) and (0, n) is pinned just below the 
Fermi energy. [I] |3j This phenomenon is unexpected in 
the band calculations and it is considered to be closely 
related to the opening of the pseudogap on the Fermi 
surface (FS) which is also an extraordinary feature 

in high-T c cuprates. This anomalous nature of the FS 
will be the direct evidence for the non-Fermi liquid be- 
havior. It is thus an interesting issue to study the FS in 
the presence of strong correlation. 

The effect of the flat band and the geometry of the 
FS can be taken into account by using the t-t'-J model 
or the t-t'-U Hubbard model in which the next-nearest 
neighbor hopping term t' is introduced as a fitting pa- 
rameter. If one chooses t' < 0, the FS centered at (tt, 
7t) observed experimentally j|H can be reproduced in 
the tight-binding model. However high temperature ex- 
pansion studies on the momentum distribution function 
for the t-J model [|| have shown that the FS is similar 
to that with t' < even though the i'-term is absent in 
the Hamiltonian. On the other hand, the conventional 
mean-field theories, such as slave-boson theory, simply 
give the FS with t' = 0. Therefore the strong correlation 
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which is not included in the mean-field theories will be 
the origin of the change of the FS geometry. 

Here we study this problem from a different point of 
view. Since the calculation in Ref. || is carried out in 
the high temperature region, it is not clear whether or 
not the FS deforms down to zero temperature. To study 
the FS of the ground state is generally very difficult. 
The exact diagonalization study of small clusters does 
not give enough resolution in the k space. The quan- 
tum Monte Carlo simulations have been often useless for 
the two-dimensional t-J model because of the minus sign 
problem. Therefore we use the variational Monte Carlo 
(VMC) method in this paper, which is free from the limi- 
tation of the system size as well as from the sign problem. 
The VMC method treats exactly the constraints of no 
doubly occupied sites and gives accurate estimates of the 
expectation values such as the variational energies and 
the momentum distribution functions. 

Although it is a variational theory, the VMC method 
is powerful to see whether some kind of symmetry 
breaking takes place or not. In this paper we ex- 
amine the Gutzwiller-projected d-wave superconducting 
state which contains an additional variational parameter 
t' v corresponding to the next-nearest neighbor hopping 
term. We can safely discuss the relative energy differ- 
ence between the variational states with and without t' v , 
although the absolute values of the variational energies 
can still be lowered. We find that the wave function with 
t' v ~ —0.1 has the lowest variational energy even though 
the Hamiltonian does not contain i'-term. This means 
that the FS deforms spontaneously. The momentum dis- 
tribution function n(k) calculated in the optimized wave 
function is consistent with that in the high temperature 
expansion. Our method gives an independent and com- 
plementary support of the result that the deformation 
of the FS is a distinctive feature of strongly correlated 
electron systems. 

In addition to this, we can identify the physical origin 
of the FS deformation in our variational approach. We 
show the relation between the energy gain and the van 
Hove singularity. It has been argued that a remarkable 
enhancement of SC correlation is achieved if the van Hove 
singularity is close to the Fermi energy. |n^,[llj Our re- 
sults show some similarity to this picture. Furthermore, 
by comparing the obtained results with the Gutzwiller 
approximation, we can see that the finite t' v is caused 
solely by the Gutzwiller projection. 

We use the two-dimensional t-J model on a square lat- 
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tice, 



H=-tJ2 P G(4a^ + h.C.)P G + JJ2^ 1 - S j} (1) 

(ij)<r {i'j) 

where (ij) represents the sum over the nearest-neighbor 
sites. c ia (ciu) is a creation (annihilation) operator of a 
(t or I) electron at i-site and Si = ^(str)^^. The 
Gutzwiller's projection operator Pq is defined as Pq = 
ILi(l — niihii), which prohibits the doubly occupied sites. 
We set J/t = 0.3. 

We use a Gutzwiller-projected mean-field type wave 
function Pg-FWJ^o) as a trial state with fixing the num- 
ber of electrons N c through Pjv e - The state is written 
as 

PgPn.\4>o) = PgPn Y[{u k + v k cl r ct kl )\0) 
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|0), (2) 



where v k /u k = A k /(e k - M + v(^^F+^f) and a «i 
is a Fourier transform of v k / u k . 

Usually e k is chosen to be e k = — 2(cosfc x + cosfc y ), 
which is in accordance with the Hamiltonian ([!]). How- 
ever in this paper, we introduce an additional variational 
parameter t' v which changes the FS of the variational 
state. We assume e k and as 

e k = —2 (cos k x + cos k y ) — 4i v cos k x cos k v , 
Afe = 2Ad(cosk x — cosky). (3) 

The present wave function contains three variational pa- 
rameters t' v , \x and Ad- 

At first, using the above wave function we calculate 
the variational energy E var of the Hamiltonian (^) 



^Pn^PqHPgPnM 

bo\PN P G P G PNM ' 



(4) 



by means of the VMC method. The distribution of the 
wave vectors k is determined in the periodic boundary 
conditions in the x direction and in the antiperiodic ones 
in the y direction so as to avoid the gap node of d-wave 
superconductivity. Although the results in the 10x10 
square lattice are mainly shown in the following, we also 
calculate larger sizes up to 20 x 20 to investigate the size 
dependence. 
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FIG. 1. Ad dependence of E val /t for t' v = (open squares), 
—0.05 (filled squares), —0.1 (open circles) and —0.15 (filled 
circles). Each point has been minimized with respect to /i. 
System size is 10x10 and the number of Monte Carlo samples 
is 2.5 x 10 5 . 

Figure [l| shows the dependence of E val - for various 
values of t' v at the doping rate S = 0.12. Apparently 
t' v ~ —0.1 gives the lowest variational energy. Since the 
Hamiltonian does not contain next-nearest neighbor hop- 
ping terms, the present result means that the shape of 
the FS of the ground state is different from that of the 
non-interacting Hamiltonian. We have also checked that, 
if the Hamiltonian has the next-nearest neighbor hopping 
term t', the optimized variational state has t' v , whose am- 
plitude is larger than t' . 

The most significant effect of this result appears in 
the shape of momentum distribution functions. Figure ^ 
shows a contour map of the gradient of the momentum 
distribution function |Vfen(fe)| for t' v = —0.1 calculated 
on the 20 x 20 square lattice. Although we have used the 
optimized variational parameters A^ and /i on the 10 x 10 
lattice, it is justified because their size dependeces are 
negligible. Brighter areas in Fig. correspond to the 
momentum k with larger values of |Vfcn(fe)|. Although 
we cannot specify exactly the location of the FS due to 
the d-w&ve SC gap, we expect that the FS lies close to 
the area where |Vfcn(fe)| is large. 

Our result of momentum distribution function is sim- 
ilar to that obtained in high temperature expansion by 
Putikka et al. || Since we take an opposite approach to 
high temperature studies, i.e., in the zero temperature, 
it is confirmed that the FS shown in Fig. || is an intrinsic 
feature of the t-J model. Note here that the smearing 
of the FS around (n, 0) in our calculation is due to the 
d-w&ve SC gap. This suggests that the similar smearing 
observed in Ref. || may be due to the pseudogap with d- 
wave symmetry, in addition to the smearing due to finite 
temperature. 

For the wide range of doping S — 0.04 ~ 0.20, we 
find that E var is minimized around t' v ~ —0.1 and the 
chemical potential /i ~ —0.5 ± 0.05. Because of the in- 
sensitiveness of fi as a function of doping, the area of the 
momentum space enclosed by the FS is also insensitive 
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to the doping rate. This result supports the violation 
of the Luttinger theorem suggested by Putikka et al. |jj 
Actually n(k) at the doping S = 0.2 in Ref. M is very 
close to our results in Fig. ^ at S = 0.12. 
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FIG. 2. Full Brillouin zone plot of |Vfen(fc)| for the doping 
rate 5 = 0.12 and J ft = 0.3. The areas with large values of 
|Vfcri(fc)| are highlighted. Note that the FS is deformed even 
though the Hamiltonian does not contain the next-nearest 
neighbor hopping terms. The smearing near (tt, 0) and (0, tt) 
is due to the presence of the d-wave SC order parameters. 

Figure || shows the energy difference between the value 
at t' v = and at t' v = —0.1, i.e. the energy gain due to the 
finite t' v , for several system sizes. Although the Monte 
Carlo results scatter a little, there is apparently a ten- 
dency that the energy gain due to the finite t' v becomes 
maximum around 5 = 0.12. As we increase the system 
size, the energy gain slightly decreases, but it will remain 
finite in the thermodynamic limit. 
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FIG. 3. Doping dependence of the energy gain at t' v — —0.1 
compared with t' v = 0. The system sizes are 10 x 10 (open 
squares), 12 x 12 (filled squares) and 14 x 14 (open circles). 
The number of samples is more than 10 J . 

Let us discuss here the relation between the energy 
gain and the van Hove singularity. For the optimized 
value t' v ~ —0.1, ek at k = (tt, 0) becomes —0.4. On the 
other hand, the optimized chemical potential fi is around 
—0.5. This means that the suddle point or the position of 
the flat band near k — (tt, 0) is very close to the chemical 
potential. Since the d-wave SC gap has a maximum at 
(tt, 0), the enhancement of the density of states near the 
Fermi energy due to the van Hove singularity causes the 
energy gain. Actually, if we assume = 0, the lowest 
energy is achieved at t' v = 0. This indicates that the FS 
deforms itself so as to fix the van Hove singularity to the 
chemical potential in the presence of the d-wave SC gap. 
This looks consistent with the mechanism of SC due to 
the van Hove singularity. 
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If we use a Hamiltonian H t -j without projection oper- 
ator and the mean-field wave function \(j)o), it is apparent 
that the variational energy {^o\Ht-j\4>o) / (^>o\^o) is min- 
imized at t' v — 0. Therefore the energy gain due to the 
non-zero value of t' y is solely from the Gutzwiller's pro- 
jection operator. In order to clarify the effect of the pro- 
jection, we examine the Gutzwiller approximation [|12| , 
in which the effect of constraints are taken into account 
by statistical weighting factors. For the t-J model, we 
have 



(^oIPg^gI^o) 

and 

{MPgPg\^) 



(5) 



= 9s(4>o\Si ■ Sj\(j> ) = g s {Si ■ Sj) , 

(G) 



where gt and g s are the renormalization factors due to the 
projection. In the simplest Gutzwiller approximation, 
g t and g s are constant, i.e., g t = 25/(1 + 6) and g s — 
4/ (l+(5) 2 . In this case, the Gutzwiller projection does 
not alter the mean-field results. However it was recently 
shown that the dependence of the renormalization factors 
on the expectation values, such as x — { c \ a c j<?)o &ud A = 
( if jO ' pl & y s a crucial role in evaluating the variational 
energies. [p^|-^6[ If we use this Gutzwiller approximation, 
we can show that 



SE V 



d(H) 
dx' 



St' + SNt'Jt' 



(7) 



where x' — ( c L c jct)o with (ij) being the next-nearest 
neighbor sites and 
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d(H) _ dg, 



(yV to) 

(8) 

The first term on the r.h.s. of eq. (0) is linear with respect 
to t' v so that -Evar is minimized at a finite value of t' v which 
satisfies 



t' 



1 d(H) 
'8N d X ' ' 



(9) 



Apparently the renormalization factors gt and g s due to 
the projection operator and their nonlinear dependence 
on x' are the origin of the spontaneous deformation of 
the FS. These phenomena cannot be found in the mean- 
field theories. The explicit calculations will be published 
elsewhere. 

In summary, we investigated the shape of the FS in the 
two-dimensional t-J model by means of the VMC calcu- 
lation introducing an additional variational parameter t' v . 
We found that the variational energy is minimized around 
t' v ~ —0.1 for various doping rates. The system size de- 
pendence indicate this effect is realized even in the ther- 
modynamic limit. The magnitude of the energy gain is 
large enough compared with other VMC studies. For ex- 
ample, the energy difference between the pure d-wave SC 
phase and the coexistent phase of AF and d-wave SC is 
comparable to the present energy gain at the doping rate 
6 - 0.08. f5|Jl7| Then we have clarified the origin of the 
energy gain by examining the van Hove singularity and 
the effect of the projection using the Gutzwiller approx- 
imation. Combining our results at zero temperature and 
those in high temperature expansion, we consider that 
the FS deformation is the most significant phenomenon 
in the presence of strong correlation. 

Since the nesting property of the FS becomes worse in 
the present wave function than in the original t-J model, 
the coexistence of AF and d-wave SC near half-filling 
Ji5| , ft7t will be suppressed. Alternatively we expect some 
incommensurate AF correlations which were unexpected 
in the t-J model. This is presumably related to the stripe 
state observed experimentally. 
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